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Within a nonlinear chiral Lagrangian framework the correlations between the quark and glueball 
\^ • admixtures of the isosinglet scalar mesons below 2 GeV and the current large uncertainties on the 

mass of the /o(600) and the /o(1370) are studied. The framework is formulated in terms of two 
scalar meson nonets (a two-quark nonet and a four-quark nonet) together with a scalar glueball. It 
£Sj ■ is shown that while some properties of these states are sensitive to the mass of /o(600) and /o(1370), 

I ' several relatively robust conclusions can be made: The /o(600), the /o(980), and the /o(1370) are 

admixtures of two and four quark components, with /o(600) being dominantly a non-strange four- 
quark state, and /o(980) and /o(1370) having a dominant two-quark component. Similarly, the 
f***) ' /o(1500) and the /o(1710) have considerable two and four quark admixtures, but in addition have 

£\J . a large glueball component. For each state, a detailed analysis providing the numerical estimates 

of all components is given. It is also shown that this framework clearly favors the experimental 
£N| ' values: m Gxp - [/o(600)] < 700 MeV and m oxp - [/ (1370)] = 1300-1450 MeV. Moreover, an overall fit 

^ ' to the available data shows a reciprocal substructure for the /o(600) and the /o(1370), and a linear 

, correlation between their masses of the form m[/o(1370)] = 0.29 m[fo (600)] + 1.22 GeV. The scalar 

■ glueball mass of 1.5—1.7 GeV is found in this analysis. 

T— I ' 

t*"- ' PACS numbers: 13.75.Lb, 11.15.Pg, 11.80.Et, 12.39.Fe 
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I. INTRODUCTION 

Exploring the properties of the scalar mesons is known to be a non trivial task in low energy QCD. This is due to 
O .' issues such as their large decay widths and overlap with the background, as well as having several decay channels over 
a tight energy range [1] . Particularly, the case of I = states is even more complex due to their various mixings, for 
' example, among two and four quark states and scalar glueballs. Below 1 GeV, the well-established experimental states 
are [1]: the / (980) [I = 0] and the ao(980) [I = 1], together with states that have uncertain properties: the /o(600) or 
a [I = 0] with a mass of 400-1200 MeV, and a decay width of 600-1000 MeV, and the ^(800) or k [I = 1/2] which 
is not yet listed but cautiously discussed in PDG. The existence of the Xq(800) has been confirmed in some theoretical 
models, while has been disputed in some other approaches. In the range of 1—2 GeV, the listed scalar states are [1]: 
K*(U30) [I = 1/2]; a (1450) [I = 1]; and / (1370), /o(1500), / (1710) [I = 0]. The isodoublet and isotriplet states, 
are generally believed to be closer to qq objects, even though some of their properties significantly deviate from such 
a description. Among the heavier isosinglet states, the /o(1370) has the largest experimental uncertainty [1] on its 
mass (1200—1500 MeV) and decay width (200—500 MeV), and other states in the energy range of around 1.5 GeV 
(or above) seem to contain a large glue component and maybe good candidates for the lowest scalar glueball state. 

A simple quark- antiquark description is known to fail for the lowest-lying scalar states, and that has made them 
the focus of intense investigation for a long time. Several foundational scenarios for their substructures have been 
considered, including MIT bag model [2], KK molecule [3] and unitarized quark model [4], and many theoretical 
frameworks for the properties of the scalars have been developed, including, among others, chiral Lagrangian of refs. 
[5-17], upon which the present investigation is based. Many recent works [18-29] have investigated different aspects 
of the scalar mesons, particularly, their family connections and possible description in terms of meson nonet (s), which 
proved a comparison with this work. 

Various ways of grouping the scalars together have been considered in the literature. For example, in some ap- 
proaches, the properties of the scalars above 1 GeV (independent of the states below 1 GeV) are investigated, whereas 
other works have only focused on the states below 1 GeV. There are also works that have investigated several possi- 
bilities for grouping together some of the states above 1 GeV with those below 1 GeV. 
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In addition to various ways that the physical states may be grouped together, different bases out of quark-antiquark, 
four quark, and glueball have been considered for their internal structure. For example, in a number of investigations 
the 1 = states above 1 GeV are studied within a framework which incorporates quark-antiquark and glueball 
components. However, in lack of a complete framework for understanding the properties of the scalar mesons it seems 
more objective to develop general frameworks in which, a priori, no specific substructure for the scalars is assumed, 
and instead, all possible components (quark-antiquark, four-quark, and glueball) are considered. The framework of 
the present work considers all such components for the 1 = states, and studies the five listed 1 = states below 2 
GeV [1]. 

Besides the generality of the framework, there are supportive indications that the lowest and the next-to- lowest lying 
scalar states have admixtures of quark-antiquark and four-quark components. For example, it is shown in [10] that the 
ao(1450) and the -ftTg(1430) have a considerable two and four quark admixtures, which provides a basis for explaining 
the mass spectrum and the partial decay widths of the I = 1/2 and 1=1 scalars [Kq (800), a (980), Kq (1430) and 
ao(1450)]. It then raises the question that if the 1=1/2 and 1=1 scalars below 2 GeV are admixtures of two and 
four quark components, why should not the 1 = states be blurred with such a mixing complexity? Answering this 
question is the main motivation of this paper, and we will see that this framework shows that there is a substantial 
admixture of the two and the four quark components for the 1 = scalars below 2 GeV. 

Motivated by the importance of the mixing for the 1=1/2 and / = 1 in this framework, the case of the 1 = states 
was initially studied in [15, 16] in which mixing with a scalar glueball is also included. In [15] the parameter space 
of the 1 = Lagrangian, which is already constrained by the properties of the 1=1/2 and 1=1 states in ref. [10], 
is studied using the mass spectrum and several two-pseudoscalar decay widths and decay ratios of the 1 = states. 
The present work extends the work of [15] by investigating in detail the effect of the mass uncertainties of the /o(600) 
and /o(1370) on the components of the 1 = states below 2 GeV. It also provides an insight into the likelihood of 
the mass of / (600) and / (1370) within their experimental values [1] in ranges 400-1200 MeV and 1200-1500 MeV, 
respectively. In addition a linear correlation between these two masses is predicted by the model. 

Specifically, we will numerically analyze the mass spectrum and perform an inverse problem: knowing the mass 
of the physical states (or in the case of the / (600) and the /o(1370) a wide experimental range for their masses) 
we search for the parameter space of the Lagrangian (which is formulated in terms of a two quark nonet and a four 
quark nonet) and find solutions that can reproduce these masses. Once the parameter space is determined it provides 
information about the properties of the two nonets. The results are summarized in Fig.l. We will show how the 
present model describes the 1 = scalar states in terms of a four-quark nonet N which lies in the range of 0.83-1.24 
GeV, together with a two-quark nonet N' in the 1.24-1.38 GeV range, and a scalar glueball that this model predicts 
in the range of 1.5-1.7 GeV (figure 1). The mass range of the two quark scalar nonet N' is qualitatively consistent 
with the expected range of 1.2 GeV from spectroscopy of p-wave mesons. We will see that this analysis shows that, 
similar to the case of I = 1/2 and 1=1 scalar mesons, the 1 = scalars have a significant admixtures of two and four 
quark components, and in addition the /o(1500) and the /o(1710) have a dominant glueball content. The dominant 
component (s) of each state is summarized in Fig. 1. 

The underlying framework of the present work has been previously applied in analyzing numerous low-energy 
processes that involve scalar mesons, and consistent pictures have emerged. The existence of the /o(600) (or a) and 
the Kq(800) (or k) and their properties have been investigated in refs. [5, 6]. The lowest-lying nonet of scalars and a 
four-quark interpretation of these states is studied in [7] and applied to irr/ scattering in [9] and several decays such 
as rj' — > 777T7T [8] and r/ => Sir [14] and radiative <j> decays [12]. In addition to ref. [10], the mixing between a two quark 
and a four quark nonets is also investigated in [11] and [17]. 

We describe the theoretical framework in Sec. 1, followed by the numerical results in Sec. 2, and a short summary 
in Sec. 3. 



II. THEORETICAL FRAMEWORK 
A. Mixing Mechanism for 1=1/2 and 1=1 Scalar States 

In ref. [10] the properties of the I = 1/2 and 1=1 scalar mesons, k(900), ifg(1430), a (980) and a (1450), in a 
nonlinear chiral Lagrangian framework is studied in detail. In this approach, a qqqq nonet N mixes with a qq nonet N' 
and provides a description of the mass spectrum and decay widths of these scalars. The Kq(1A30) and the ao(1450), 
are generally believed to be good candidates for a qq nonet [1] , but some of their properties do not quite follow this 
scenario. For example, in a qq nonet, isotriplet is expected to be lighter than the isodoublet, but for these two states 
[1]: 



m [a (1450)] = 1474 ± 19 MeV > m [1^(1430)] = 1412 ± 6 MeV 
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FIG. 1: Prediction of the present model for the formation of the 7 = scalar states below 2 GeV. On the right, the physical 
states are shown, with the solid lines representing the experimentally well established masses, and the height of the two 
boxes representing the prediction of the present model for the uncertainties of the mass of /o(600) and /o(1370), which are 
approximately in ranges 0.4 - 0.7 GeV and 1.3 - 1.45 GeV, respectively. On the left, the "bare" states are shown, with the 
dashed lines representing the four-quark nonet N, the solid lines representing the two-quark nonet N' and the box representing 
the scalar glueball. The mass of the bare states in the two nonets (from bottom to top) are: m{udud) — 0.83 GeV, m(dsud) = 
1.06 GeV, m(sdds + suus)/V2) = m(uu + dd)/V2) = 1.24 GeV, m(us) = 1.31 GeV and m(ss) = 1.38 GeV. The uncertainty 
of the glueball mass (shown by the height of the box, approximately between 1.5 GeV to 1.7 GeV) is due to the uncertainty of 
the mass of /o(600) and /o(1370). The arrows show the dominant component(s) of each physical state. More details are given 
in Figs. 7 and 8. 



Also their decay ratios given in PDG [1] do not follow a pattern expected from an SU(3) symmetry (given in paren- 
thesis): 

- 0.92 ±0.12 (1.51) 



r[# *- 


-» ttK] 


r[a - 


KK] 


r [a - 


-» 7T77] 


r [a - 


-> Trr/] 


r> 


-» TT?7] 



15 ±0.16 (0.16) (2) 

These properties of the Kg (1430) an d the ao(1450) are naturally explained by the mixing mechanism of ref. [10]. The 
general mass terms for the 7=1/2 and the I — 1 states can be written as 

^maJs 2 ' 1 = -aTr(NN) - bTr(NNM) - a'Tr(N'N') - b'Tr(N'N'M) (3) 

where A4 = diag(l, l,x) with x being the ratio of the strange to non-strange quark masses, and a,b,a' and b' arc 
unknown parameters fixed by the unmixed or "bare" masses (denoted below by subscript "0"): 



m 2 [a ] = 2{a + b) m 2 [a'^ = 2(a' + b') 

m 2 [K ] = 2a+(l + x)b m 2 [K' ] = 2a' + (l + x)b' 



2r^„i — o„ _l n _l <^/, ™2rV"i _ nr.i _i_ ft _i_ „\h' (4) 



As TV is a four-quark nonet and N' a two-quark nonet, we expect: 

to 2 [K ] < m 2 [ao] < m 2 [a' ] < m 2 [K' ] (5) 

In fact, this is how we tag N and N' to a four-quark and a two-quark nonet, respectively. Introducing a general 
mixing 

<^ Al = -7Tr(AW') (6) 
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it is shown in [10] that for 0.51 < 7 < 0.62 GeV 2 , it is possible to recover the physical masses such that the 
"bare" masses have the expected ordering of (5). In this mechanism, the "bare" isotriplet states split more than 
the isodoublcts, and consequently, the physical isovector state ao(1450) becomes heavier than the isodoublet state 
.KTq(1430) in agreement with the observed experimental values in (1). The light isovector and isodoublet states 
are the a (980) and the X o *(800). With the physical masses m[a (980)] = 0.9835 GeV, m[K$ (800)] = 0.875 GeV, 
m[ao(1450)] = 1.455 GeV and m[-Ko(1430)] = 1.435 GeV, the best values of mixing parameter 7 and the "bare" 
masses are found in [10] 

m ao = m a < o = 1.24 GeV, m Ka = 1.06 GeV, m K > = 1.31 GeV, 7 = 0.58 GeV 2 (7) 

These parameters are then used to study the decay widths of the 7 = 1/2 and 1=1 states [10]. A general Lagrangian 
describing the coupling of the two nonets N and N' to two-pseudoscalar particles are introduced and the unknown 
Lagrangian parameters are found by fits to various decay widths. 

In summary, the work of ref. [10] clearly shows that properties of the lowest and the next-to-lowest 1=1/2 and 
1 = 1 scalar states can be described by a mixing between a quark-antiquark nonet and a four quark nonet. It is 
concluded that the 1=1 states are close to maximal mixing (i.e. ao(980) and do(1450) are approximately made of 
50% quark-antiquark and 50% four-quark components), and the 1=1/2 states have a similar structure with .Kg (800) 
made of approximately 74% of four-quark and 26% quark-antiquark, and the reverse structure for the .Kg (1430). 
Now if the lowest and the next-to-lowest 1=1/2 and 1=1 states have a substantial mixing of quark-antiquark 
and four-quark components, then it seems necessary to investigate a similar scenario for the 1 = scalars, which, in 
addition, can have a glueball component as well. In other words, if the lowest and the next-to-lowest scalars are going 
to be grouped together, then all components of quark-antiquark, four-quark and glueball for the 1 = states should 
be taken into account. The case of I = states will be discussed in next section. 



B. Isosinglet States 



The 1 = scalars have been investigated in many recent works, including the general approach with two and four 
quark components as well as a glueball component of refs. [15, 16]. Within the framework of ref. [10], an initial 
analysis is given in [15] in which the mass and the decay Lagrangian for the 1 = scalar mesons are studied in some 
details. The general mass terms for nonets N and N', and a scalar glueball G can be written as: 

= ^aJs' 1 ~ cTt(N)Tt(N) - dTr(N)Tr(NM) 

-cTr(N')Ti(N') - d'Ti(N')Tr(N'M) - gG 2 (8) 

The unknown parameters c and d induce "internal" mixing between the two 1 = flavor combinations [(TV-j 1 + 7Vf ) / \/2 
and iVg] of nonet N. Similarly, d and d 1 play the same role in nonet N' . Parameters c, d, c' and d' do not contribute 
to the mass spectrum of the 1=1/2 and 1=1 states. The last term represents the glueball mass term. The term 

Zmasi 2 ' 1 is imported from Eq. (3) together with its parameters from Eq. (7). 

The mixing between N and N', and the mixing of these two nonets with the scalar glueball G can be written as 



C 1 





mix 



c 



1=1/2,1 



pTr(N)Tr(N') - eGTr (N) - /GTr (N') 



where the first term is given in (6) with 7 from (7). The second term does not contribute to the I 
and in special limit of p — > —7: 

- 7 Tr(AW') - pTr(N)Tr(N') = ~fe abc e ade N*N' e c 



(9) 

1/2, 1 mixing, 
(10) 



which is more consistent with the OZI rule than the individual 7 and p terms and is studied in [21]. Here we do not 
restrict the mixing to this particular combination, and instead, take p as a a priori unknown free parameter. Terms 
with unknown couplings e and / describe mixing with the scalar glueball G. As a result, the five isosinglets below 2 
GeV, become a mixture of five different flavor combinations, and their masses can be organized as 



mass ' ^rnirr 



iF M 2 F = -iFML s F 



(11) 



with 



F n = 



N! 



\ 



N 



,3 



V 



(N'i + N'l)/V2 



G 



( udud \ 

(sdds + suus) I v2 
ss 

(uu + dd)/V2 
\ G J 



ffo NS 


\ 


fo S 




f/S 
J 




f/NS 
J 




\G 


J 



(12) 



■5 



where the superscript NS and S respectively represent the non-strange and strange combinations. F contains the 
physical fields 



/ ^(550) \ 
/ (980) 
/o(1370) 
/o(1500) 

V /o(1710) / 



where K 1 is the transformation matrix. The mass squared matrix is 
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(14) 



in which the value of the unmixed / = 1/2, 1 masses, and the mixing parameter 7 are substituted in from (7). We see 
that there are eight unknown parameters in (14) which are c, d, c', rf', g, p, e and /. We will use numerical analysis 
to search this eight dimensional parameter space for the best values that give closest agreement with experimentally 
known masses. Particularly, we will study in detail the effect of the large uncertainties on the mass of the /o(600) 
and /o(1370) on the resulting parameters. 



III. MATCHING THE THEORETICAL PREDICTION TO EXPERIMENTAL DATA 

To determine the eight unknown Lagrangian parameters (c, c', d, d' , e, /, g and p), we input the experimental masses 
of the scalar states. Out of the five isosignlet states, three have a well established experimental mass [1]: 

m cx P[/ (980)] = 980±10MeV 
m cxp [/ (1500)] = 1507±5MeV (15) 
m cx P[/o(1710)] = 1713±6MeV 

However, the experimental mass of the /o(600) and the /o(1370) have very large uncertainties [1]: 

TO cxp -[/ (600)] = 400 -» 1200 MeV 
m oxp -[/o(1370)] = 1200 ^ 1500 MeV (16) 

We search for the eight Lagrangian parameters (which in turn determine the mass matrix (14)) such that three of 
the resulting eigenvalues match their central experimental values in (15), and the other two masses fall somewhere in 
the experimental ranges in (16). This means that out of the five eigenvalues of (14), three (eigenvalues 2, 4 and 5) 
should match the fixed target values in (15), but two (eigenvalues 1 and 3) have no fixed target values and instead 
can be match to any values in (16). Therefore, to include all possibilities for the experimental ranges in (16), we 
numerically scan the m exp ' [/o (600)]— m cxp [/o (1370)] plane over the allowed ranges, and at each point we fit for the 
eight Lagrangian parameters such that the three theoretically calculated eigenvalues (2, 4 and 5) match the fixed 
target masses in (15), and the other two eigenvalues (1 and 3) match the variable target masses at the chosen point 
in this plane. At a given point in the m cxp - [/ (600)]— m cxp [/ (1370)] plane, we measure the goodness of the fit by 
the smallness of the quantity 

| m thoo. _ CXP-I 

x(m oxp -[/ (600)], m cxp -[/ (1370)]) = £ U (17) 

1 l 

where i = 1 • • -5 correspond to the five / — states in ascending order (for example m c * v ' — m cxp - [/ (600)], • • •). The 
value of xx 100 gives the overall percent difference between theory and experiment. The total of 14,641 eight-parameter 
fits were performed over the target points in the m exp ' [/o (600)]— m cxp [/o (1370)] plane (the overall completion of the 
numerical analysis of this article required several months of computation time on a XEON dual-processor workstation). 
This procedure creates a three dimensional graph of x as a function of m cxp [/ (600)] and m exp - [/ (1370)]. The overall 
result of the fits for the experimentally allowed range of masses are given in Fig. 2, in which the projection of \ onto 
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FIG. 2: Projection of x onto the \ ~ m cxp ' [/ (600)] plane, and onto the x ~ m cxp [/ (1370)] plane. The left figure shows that 
values of m exp ' [/o(600)] > 700 MeV result in a significant disagreement between theory and experiment. The right figure shows 
that m cxp '[/ (1370)] outside of the range 1300-1450 MeV is not favored by the model. 



the x~ m exp [/o(600)] plane and onto the % — m exp [/o(1370)] plane are shown. We can easily see that the experimental 
mass of /o(600) above 700 MeV, and the experimental mass of /o(1370) outside the range of 1300 to 1450 MeV are 
not favored by this model. 

More refined numerical analysis narrows down the favored regions to the ranges shown in Fig. 3. We see that the 
lowest value of \ occurs within the ranges: 

TO oxp -[/ (600)] = 500 -^600 MeV 
m cxp [/ (1370)] = 1350 ^ 1400 MeV (18) 

In this region, \ exhibits a series of local minima, also shown in a more refined numerical work in Fig. 4. Although 
X has its lowest values in the regions given in (18) we see that it only increases by approximately 0.01 outside of this 
region and can be comparable to the theoretical uncertainties of this framework. Therefore, for a more conservative 
estimate we take into account all local minima with x < 0.02. This confines the experimental masses to the ranges: 

m cxp [/o(600)] = 400 700 MeV 
m oxp [/o(1370)] = 1300 ^ 1450 MeV (19) 

The local minima of \ in this range are shown in Fig. 5 in which the projection of x onto the 
TO oxp [/o(600)]— m cxp [/o(1370)] plane is given. In the gray region, there is an overall disagreement of less than 
5% between theory and experiment. The local minima (with \ < 0.02) are shown with dots, together with a linear 
fit that shows the correlation between the mass of the /o(600) and the /o(1370): 

m[/ (1370)] = 0.29 m[/ (600)] + 1.22 GeV (20) 

At all local minima (with x < 0.02), a detailed numerical analysis is performed and the eight Lagrangian parameters 
are determined. The large uncertainties on m[/o(600)] and m[/o(1370)] do not allow an accurate determination of these 
parameters. In this work we study the correlation between these uncertainties and determination of the Lagrangian 
parameters. For orientation, let us first begin by investigating the central values of 



m oxp [/ (600)] = 558 MeV 
m exp [/o(1370)] = 1373 MeV 



(21) 
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FIG. 3: Projection of \ ont o the \ — m cxp ' [/o(600)] plane, and onto the \ ~ m° xp ' [/o(1370)] plane. The left figure shows a 
series of local minima in x with the lowest values in the range 500 MeV < m exp ' [/o(600)] < 600 MeV. The right figure shows 
that the present model predicts that m exp ' [/o(1370)] is confined within a smaller range of 1300 MeV to 1450 MeV. 



It is important to also notice that the central value of to[/o(600)] = 558 MeV is exactly what was first found in [5] in 
applying the nonlinear chiral Lagrangian of this work to irir scattering. At the particular point of (21), the result of 
the fit is given in Table I, and is consistent with the initial investigation of this model in rcf. [15] in which the effect 
of the mixing parameter p is only studied at several discrete points. The result given in Table I supplements the work 
of [15] by treating p as a general free parameter. The rotation matrix is 



and in turn determines the quark and glueball components of each physical state as presented in Fig. 6. Several 
general observations can be made in Fig. 6. Clearly, the /o(600) is dominantly a non-strange four-quark state with 
a substantial ss component. On the other hand the dominant structure of the / (1370) seems to be the reverse of 
the /o(600). The /o(980) seems to have a dominant non-strange two-quark component with a significant strange 
four-quark component. The /o(1500) appears to have a dominant glueball component with some minor two and four 
quark admixtures. The dominant component of the / (1710) seems to be a strange four-quark state followed by a 
glueball component. 

The scalar glueball mass can be calculated from the fit in Table I: 



which is in a range expected from Lattice QCD [30]. 

Next we will examine the effect of deviation from the central values of to[/o(600)] and m[/o(1370)] on the predictions 
given in Fig. 6. We will find that the average properties of the /o(600), /o(980) and /o(1370) remain close to those 
given in Fig. 6, but some of the properties of the /o(1500) and /o(1710) are sensitive to such deviations. 

To investigate the effect of the mass uncertainties of the /o(600) and /o(1370), we perform 8-parameter fits at each 
of the local minima of Fig. 5, and determine the variation of the Lagrangian parameters around the values in Tabic 
I. The results are summarized in Table II. We see that practically the Lagrangian parameters are quite sensitive to 



K- 1 



- 0.812 -0.065 -0.542 -0.183 0.101 " 

0.346 -0.447 0.294 0.764 -0.106 

0.418 0.198 0.768 -0.328 0.299 

-0.106 0.292 -0.174 0.403 0.844 

0.190 0.820 -0.006 0.338 -0.422 



(22) 




(23) 
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FIG. 4: Projection of x onto the x — rn c * p ' [/o(600)] plane, and onto the x — m cxp '[/o(1370)] plane, x exhibits a series of local 
minima. 
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FIG. 5: Projection of x onto the m exp ' [fo (600)]— m cxp ' [fo (1370)] plane. In the gray region there is an overall disagreement of 
less than 5% between theory and experiment. The circles represent points on this plane at which x has a local minimum and 
X < 0.02. The solid line captures the trend of these local minima with equation: m[/o(1370)] = 0.29 m[fo (600)] + 1.22 GeV. 

the mass of the /o(600) and /o(1370). The only exception is parameter g which determines the glueball mass: 

m G = y/Yg = 1.5 -» 1.7 GeV (24) 

At the local minima of Fig. 5, we compute the rotation matrix K^ 1 (defined in Eq. (13)) which maps the quark 
and glueball bases to the physical bases. The results show that K^ 1 is less sensitive to the mass of the /o(600) and 
/o(1370). For each scalar state, the quark and glueball components versus m oxp - [/o(600)] are given in Fig. 7, and the 
averaged components together with their variations are given in Fig. 8. Although the m cxp ' [Jo (600)] is still over a 
wide range of 400 to 700 MeV, we see in Fig. 7 that some of the components of the physical state are qualitatively 
stable: The admixtures of / (600), /o(980) and /o(1370), as well as some of the components of the /o(1500) and 
/o(1710) are within a relatively small range of variation. We see that the /o(600) has a dominant udud component 
and some ss content, and has almost the reverse structure of the / (1370); the / (980) is dominantly a two quark 
non-strange state (uu + dd) with a significant strange four-quark content (sdds + suus). Although sensitive to the 
mass of the /o(600) and / (1370), we see that the /o(1500) contains a dominant glueball content with some strange 
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Lagrangian Parameters 


Fitted Values (GeV 2 ) 


c 


2.32 xKT 1 


d 


-9.11 x 10~ 3 


c' 


-4.29 x 1(T 3 


d' 


-1.25 x 10~ 2 


9 


1.16 


P 


2.87 x 1(T 2 


e 


-2.13 x KT 1 


f 


6.10 x 10~ 2 



TABLE I: Fitted values of the Lagrangian parameters for m[/ (600)] = 558 MeV and m[/ () (1370)] = 1373 MeV. 



Lagrangian Parameters 


Average (GeV 2 ) 


Variation (GeV 2 ) 


c 


1.71 x 10 _1 


(0.37 -> 2.60) x 10" 1 


d 


-7.32 x 10~ 3 


(-14.44 -> -2.63) x 10~ 3 


d 


-2.87 x 1(T 3 


(-5.32 -> -0.64) x 10~ 3 


d! 


-0.97 x 10~ 2 


(-1.41 -> -0.27) x 10~ 2 


9 


1.24 


1.15 -> 1.44 


P 


4.7 x 10~ 2 


(1.31 -> 9.14) x 10~ 2 


e 


-2.10 x 10 _1 


(-3.02 -> -1.76) x 10" 1 


f 


3.74 x 10~ 2 


(1.15 -> 11.40) x 10~ 2 



TABLE II: The effect of the uncertainties of m exp ' [/o(600)] and m exp ' [/o(1370)] on the Lagrangian parameters. The averaged 
value of the parameters (second column) is compared with their (asymmetric) range of variation (third column) . Only parameter 
g (which determines the glueball mass) is relatively insensitive to m exp ' [/o(600)] and m cxp ' [/o(1370)]. 



four-quark and non-strange two-qurak admixtures; and the /o(1710) is a dominant strange four-quark state with a 
comparable glueball component. 

IV. SUMMARY AND CONCLUSION 

We studied the 1 = scalar mesons below 2 GeV using a nonlinear chiral Lagrangian which is constrained by the 
mass and the decay properties of the / = 1/2 and 1=1 scalar meson below 2 GeV [Kq(800), Kq(1430), a (980) 
and ao(1450)]. This framework provides an efficient approach for predicting the quark and glueball content of the 
scalar mesons. The main obstacle for a complete prediction is the lack an accurate experimental input for the mass 
of /o(600) and /o(1370). Nevertheless, we showed that several relatively robust conclusions can be made. We showed 
that the /o(600), the / (980) and the /o(1370) have a substantial admixture of two and four-quark components with 
a negligible glueball component. The present model predicts that the /o(600) is dominantly a non-strange four-quark 
state; the /o(980) has a dominant non-strange two-quark component; and the /o(1370) has significant ss admixture. 
We also showed that this model predicts that the /o(1500) and the /o(1710) have a considerable two and four quark 
admixtures, together with a dominant glueball component. The current large uncertainties on the mass of /o(600) 
and /o(1370) do not allow an exact determination of the glueball components of the /o(1500) and the /o(1710), but it 
is qualitatively clear that the glueball components of these two states are quite large. In addition the present model 
predicts that the glueball mass is in the range 1.5—1.7 GeV (Eq.(24)). 

The main theoretical improvement of the model involves inclusion of higher order mixing terms among nonets N 
and N' and the scalar glueball: 

Tr (MNN') + Tr (MN'N) , Tr{MN)Tr(N') , Tr(MN')Tr{N) , e'GTr (MN) , f'GTr(MN') (25) 

However, these terms both mix the quarks and glucballs as well as break SU(3) symmetry, and therefore are of a 
more complex nature compared to the terms used in this investigation (Eqs. (6) and (9)). Investigation of such higher 
order mixing terms will be left for future works. It is also important to note that the results presented here are not 
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sensitive to the choice of 7 which determines the mixing among the 1=1/2 and 1=1 mixings in Eq. (6). This is 
shown in Fig. 9 in which the components of all / = states are plotted versus 7 2 and show that they are relatively 
stable. Also in Fig. 10 the bare masses are plotted versus 7 2 in which we see that the expected ordering (i.e. the 
lowest-lying four-quark nonet N underlies a heavier two-quark nonet N' and a glueball) which is examined in ref. 
[10] using the properties of the 1 = 1/2 and 1=1 scalar states, is not sensitive to the choice of mixing parameter 7, 
providing further support for the plausibility of the leading mixing terms considered in the present investigation. 
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FIG. 6: Percentage of the quark and glueball components of the scalar states for values of m[/o(600)] = 558 MeV and 
m[/o(1370)] = 1373 MeV. In each figure, the columns from left to right respectively represent the percentage of: udud (black), 
(sdds + suus)/\/2 (light gray), ss (dark gray), (mm + dd)/\/2 (white), and glueball (black). 
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FIG. 7: Percentage of quark and glueball components of the scalar states vs m exp " [/o(600)]: udud (filled circles), (sdds + 
suus)/\/2 (empty circles), ss (filled squares), (uu + dd)/y/2 (empty squares), and glueball (filled triangles). 
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FIG. 8: The effects of the mass uncertainties of the /o(600) and /o(1370) on the percentage of the quark and glueball components 
of the scalar states. Components 1 to 5 respectively represent udud, (sdds + suus)/V^, ss, (uu + dd)/\/2, and glueball. 
The dots represent the averaged values of each component and the error bars reflect the uncertainties of m oxp ' [/o (600)] and 
jn ° xp '[/o(1370)]. The figures show that the components of the /o(600), /o(980) and /o(1370) are not very sensitive to these 
experimental uncertainties, but some of the components of the /o(1500) and /o(1710) (such as their glueball component) are 
significantly affected by such uncertainties. 
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FIG. 9: Percentage of quark and glueball components of the scalar states vs y 2 : udud (filled circles), (sdds + suus) / \/2 (empty 
circles), ss (filled squares), (uu + dd) / y/2 (empty squares), and glueball (filled triangles). The figures show that the components 
are not very sensitive to the mixing parameter 7. 
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FIG. 10: Dependence of the bare masses on . The four-quark scalar nonet (dashed lines) lies below the quark-antiquark 
scalar nonet (solid lines) and a scalar glueball (dotted-dashed line). 



